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Abstract 

In this contribution I discuss the nuclear symmetry energy in the regime of 
hadronic degrees of freedom. The density dependence of the symmetry energy is 
important from very low densities in supernova explosions, to the structure of neutron- 
rich nuclei around saturation density, and to several times saturation density in neu¬ 
tron stars. Heavy ion collisions are the only means to study this density dependence 
in the laboratory. Numerical simulations of transport theories are used to extract the 
equation-of-state, and thus also the symmetry energy. I discuss some examples, which 
relate particularly to the high density symmetry energy, which is of particular interest 
today. I review the status and point out some open problems in the determination of 
the symmetry energy in heavy ion collisions. 

1 Introduction 

The nuclear Equation-of-State (EoS) is often taken to specify the energy of nuclear 
matter as a function of density, temperature and asymmetry. For zero temperature it 
can be written as E(p, 5) = E nm (p)-\-E sym (p)5 2 -\-..., where p is the total density of the 
system and 5 = (p n —p p ) /P the asymmetry with the neutron and proton densities. The 
term proportional to S 2 is the symmetry energy. For saturation density it is related to 
the symmetry energy term in the empirical mass formula. However, the dependence 
on density is of great importance in nuclei away from stability and in astrophysics 
in core-collapse supernovae and neutron stars, which have a large neutron excess 
and where a large range of densities from very low (in supernovae) to several times 
saturation density (in neutron stars) is involved. Predictions of mircroscopic many- 
body calculations of the symmetry energy differ widely, especially above saturation 
density pQj . Most likely the reason is the poorly known short range isovector repulsion 
[2]. Thus there are extensive efforts to determine the symmetry energy from nuclear 
structure, heavy ion collisions and astrophysical observations. 

In nuclear structure one explores the symmetry energy around saturation density 
Po, using an expansion of the form E sym = So + (L/3)(p —po)/p + ..., or a parametriza- 
tion which splits the symmetry energy into a kinetic contribution taken in the form of 
the free Fermi gas, and a potential part as a power law E sym = |ei?(p/ po) 2 ^ 3 +C(p/ po) 7 
with a parameter 7 . The slope of the symmetry energy L , or more specifically the 
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correlation between the value and the slope, So vs. L , has been extensively investi¬ 
gated using various observables, like nuclear masses, Giant and Pygmy dipole res¬ 
onances, dipole polarizabilities, neutron skin radii (difference between neutron and 
proton radii), and isobaric analog state energies [2113]. Neutron star observations, on 
the other hand, also provide access to the nuclear EoS. since it determines uniquely 
the mass-radius relation. The simultaneous measurement of these two properties is 
difficult but significant progress has been made in the last years, using sophisticated 
modelling of neutron star atmospheres and statistical analyses [5]. However, definite 
conclusions are still controversial. 

In this situation heavy ion collisions (HIC) from Fermi energies up to intermediate 
energies of several GeV per particle provide a means to investigate the nuclear EoS in 
the laboratory. In a collision nuclear matter first compresses and then expands in the 
final state, and thus different regions of density are explored. The advantage of heavy 
ion collisions is the freedom to vary the energy and the impact parameter (and thus 
the compression), and also the asymmetry of the colliding system within limits, which 
will be extended in the future with new rare isotope facilities. The difficulty is that a 
HIC is fundamentally a non-equilibrium process, and thus the detailed evolution has 
to be modelled using transport theory, on which there are still open questions. Also 
one has to identify observables which are especially sensitive to the symmetry energy 
in the presence of uncertainities of the much larger contribution of symmetric nuclear 
matter and the above-mentioned uncertainities of transport theory. 

Here I aim to give a brief overview of the present status of the investigation of 
the symmetry energy using HIC, using not only results of my collaborations but also 
those of other groups. Certainly this is not possible in any complete way in the limited 
space of this article. Thus I will select specific examples, which also point to the still 
open problems in these investigations. Recently there have been extensive reviews 
on the symmetry energy in the form of volumes collecting articles of experts M, 
review articles M, and feature articles Pll], just to mention the more recent ones, 
which contain much more information. 

2 Theoretical Considerations 

The main method to interpret HIC is a transport theory, which describes the tem¬ 
poral evolution of the one-body phase-space distribution function /( r, p, t) under the 
action of a mean field potential G(r,p), possibly momentum dependent, and 2-body 
collisions with the in-medium cross section a(f2). In a non-relativistic approach it 
reads 


dh_ 

dt 


^ + — V {r >/i 
at m 


V w c(r, p)- VW[/ ; (r, p) V' r >/ ; 
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= J2 dPjdPi'dPj'Uij^ij->i'j'(^)5(Pi + Pj - Pi' - Pj' 

x[(l - /i)(l - - fifj ( 1 - /*')(! - fj')], 


but field theoretical formulations are also widely used [§]. The indices run 

over neutrons and protons, such that these are coupled equations via the collision term 
and indirectly via the potentials. If the production of other particles is considered, like 
A’s, or Ti and K mesons, these have their own transport equations coupled through 
the corresponding inelastic cross sections. The mean field potentials U t can be derived 
from an energy functional. Mean fields and cross sections should be related through a 
theory for the in-medium effective interaction, like e.g. Brueckner theory, even though 
this is not necessarily done in many applications. The isospin effects enter via the 
differences in neutron and proton potentials and the isospin dependent cross sections, 
but they are always small relative to the dominant isoscalar effects. Thus one often 
resorts to differences or ratios of observables between isospin partners, in order to 
eliminate as much as possible the uncertainties in the isoscalar part. The decription 
of cluster production, which as seen below is often an issue in the investigations of the 
symmetry energy, in principle goes beyond the one-body description, and its proper 
treatment is still one of the questions in the application of transport theory. 

3 Isospin observables 

It may be useful to briefly summarize the different ways in the investigations of 
the symmetry energy in the different density regimes. In central reactions at Fermi 
energies densities somewhat above saturation are reached. Recently the expansion 
phase of such reactions has been studued in detail, where very low densities of about 
1/10 to 1/1000 of p 0 are attained. From isostope ratios (so-called iso-scaling) the 
symmetry energy at very low densities has been determined which is important for 
the simulation for supernova explosions. In this density regime few-body clustering 
effects become important. A theoretical investigation has shown that the symmetry 
energy, in fact, is finite a very low densities in qualitative agreement with experiment 

m- 

At Fermi energy collisions one observes the phenomenon of multifragmentation. 
The distribution of the isospin to the different fragments (’’isospin fractionation”) [8], 
and the isospin transport through the low-density neck (’’isospin diffusion”) in more 
peripheral collisions have been very useful to constrain the symmetry energy below 

do PH- 

At intermediate energies, which are the main interest here, the initial phase of the 
collision is characterized by pre-equilibrium emission of high energy particles and light 
fragments. The yield ratios of isotopic partners, like n/p or 3 He/t, contain informa¬ 
tion on the relative strengths of the neutron and proton potentials. The compression 
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Figure 1: (left) The neutron-proton ratio in 136 A"e+ 124 SVi collisions at 150 AMeV for 
different choices of the symmery energy (solid asy-soft, dashed asy-stiff) and orderings 
of the effective masses (blue, m* < m*, red m* > m*). On the right the corresponding 
ratio of tritium over ‘‘He. 

phase then determines the momentum distribution of the emitted particles, generally 
called ’’flow”, in-plane (directed) and out-of-plane (elliptic). Neutron-proton differ¬ 
ences of flow observables have been an important means to extract information on 
the high density symmetry energy. Inelastic NN collisions lead to the production of 
A resonances, which may decay into pions or lead to the production of strangeness. 

Let us briefly note that the influence of the symmetry energy has also been dis¬ 
cussed at higher energies, suggesting that the deconfinement transition may be sub¬ 
stantially influenced by the difference of the symmetry energy between the hadronic 
and partonic phases, and may, in fact, occur at lower density in asymmetric systems 

ra 


3.1 Pre-equilibrium Emission 

The neutron to proton ratio of emitted particles has first been measured at MSU 
for Sn + Sn systems at 50 AMeV, and a systematic analysis of several observables 
has yielded rather good limits on the 7 exponent around 7 & 0.6 m • In the pre¬ 
equilibrium emission at higher energies the momentum dependence of the symmetry 
potential, i.e. the proton-neutron effective mass splitting, becomes important as first 
pointed out in refs. H3JE]. Recently we have systematically studied this effect for 
nucleons and light clusters in different Xe + Sn reactions at energies between 32 and 
150 AMeV in ref.|[T5]. A result from these calculations is shown in FigJT]for central 
collisions at 150 AMeV with different stiffnesses of the symmetry energy and different 
effective mass splittings, in the left panel the n/p ratio is shown as a function of the 
transverse energy of the emitted particles, on the right the corresponding result for 
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t/ 3 He. One observes a clear pattern, namely that the stiffness of the symmetry 
energy governs the lower part of the transverse energy spectrum, where a softer 
symmetry energy yields a larger n/p ratio. The higher energy part of the spectrum 
is dominated by the effective mass ordering, where a smaller neutron effective mass 
favors the emission of neutrons and increases the ratio. A similar result has been 
obtained by Zhang et al.|16j. Thus this observable should serve as a promising probe 
to disentangle the density and momentum dependences of the symmetry potential at 
higher density. The ratio t/ 3 He in the right panel shows a very similar pattern. A 
not yet very conclusice comparison with INDRA data |T5j favors a stiff symmetry 
energy with m* > rn*. 

3.2 Flow 

The momentum distribution of the particles and fragments emitted in the final stage 
of a HIC are characterized via a Fourier series expansion of the azimuthal distribution 
as N(Q\p t) y) = A r 0 (l — vi(p t ,y)cosQ + V2(pt.,y)cos2Q + ... The first two Fourier 
coefficients, depending on the transverse momentum p t and the longitudinal rapidity 
y, are called directed and elliptic flow, respectively. Differences of flow parameters 
between isospin partners directly reflect the isospin-dependent potentials and thus 
the symmetry energy. Preliminary results from the ASYEOS experiment of the FOPI 
collaboration are shown in Fig]2]for Au+Au collisions at 400 AMeV [IT]. Data for 
the ratio of the elliptic flow of neutrons relative to hydrogen are shown together with 
calculations for a soft and a stiff symmetry energy. A best fit yields a coefficient 7 of 
about 0.76, i.e. a moderately soft symmetry energy. This is an important advance in 
trying to constrain the symmetry energy at supersaturation densities. 

3.3 Particle Production 

The n/p asymmetry of the compressed system also influences the ratio of newly 
produced particles, which thus can serve as indicators of the symmetry energy in 
the high density phase. Pions are produced predominantly via the A resonances, 
NN ^ N A and the subsequent decay A —> Nti. The ratio of the isospin partners 
tt~/ n + can thus serve as a probe of the high density symmetry energy. As analyzed 
in ref. [18] there are competing effects on the A and pion production from the isospin 
dependent mean fields and the A production threshold conditions. 

In Figj3]l have collected in the right panel results from recent theoretical analyses 
of this ratio using different models of symmetry energies, shown correspondingly in 
the right panel [19], and different program codes. They are compared to the FOPI 
data [2D] - For each model the results for two parameter sets of different stiffness are 
shown (stiffer - blue, softer - red). As is seen, the results of the different models are 
not only very different quantitatively, but even the trend with the asy-stiffness is not 
consistent. A reason may lie in different modelling of the A dynamics, and also in the 
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Figure 2: Ratio of elliptic flow of neutrons over hadrogen for forward rapidities in 
Au+Au collisions at 400 AMeV. Data of ref. mi and calculations with a symmetry 
energy characterized by exponent 7 . 


competing mean field and the threshold effects, where slightly different treatments 
might lead to large differences. This issue needs clarification in view of the sensitivity 
of the pion observables and the excellent data situation. More information should be 
gained by discussing not the energy integrated yield ratios, but the spectral behavior 
of the ratio, since different energy pions are expected to have different histories [ 21 ] . 

ft has also been suggested, that the ratio of the anti-strange kaon isospin part¬ 
ners, K°/K + could be a useful observable for the symmetry energy [22]. Indeed, 
kaon production has been one of the most useful observables to determine the EOS 
of symmetric nuclear matter. The anti-strange kaons weakly interact with nuclear 
matter and are thus a direct probe of the dense matter where they are produced. 
Theoretical analyses show similar if not larger sensitivity to the symmetry energy 
compared to pion ratios, but appropriate date do not yet exist. 

4 Discussion and Conclusions 

In Figjf] I attempt to give a summary of the present information on the density 
dependence symmetry energy E sym (p) or equivalently S(p) from HIC [3] 0]. In the 
left panel the region below saturation density po is shown. The blue hatched area is 
the result from the investigation of Sn+Sn collisions at 50 MeV from MSU [3] using 
various observables from isospin transport between different Sn isotopes. The isolated 
symbols represents information from the fits to nuclear masses or GDR energies, which 
are plotted at about 0.6po, which is an average density of nuclei where different models 
of the EoS converge. The blue-bordered areas are derived from an analysis of isobaric 
analog states which give information also on the lower densities in the surface [23] . 
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Figure 3: (right) The 7T~/tt + ratio in Au + Au collisons as a function of incident 
energy as measured by the FOPI collaboration, and as calculated by different groups, 
as indicated in the legend and discussed in the text. On the left are shown the 
corresponding models for the symmetry energy. 


When the analysis is combined with the information on the neutron skin radius of 
Pb, the constraint is still sharpened. The points in the lower left corner, labelled 
’’cluster”, come from an analysis of the very low density matter in the expansion 
phase of low energy heavy ion collisions, mentioned above m- Here the matter is 
not any more homogeneous but cluster correlations become important, which have 
the effect of making the symmetry energy finite at very low densitiy. Alltogether the 
various sources of information on the symmetry energy in this density range seem 
to converge, and they also converge with the theoretical many-body results, see e.g. 
ref. P . 

The information on the symmetry energy above saturation is shown in the right 
panel (where the low densitiy results from the right panel are shown by the red- 
bordered area), particularly the results from the neutron/hydrogen flow analysis from 
sect. 3.3, and two results from the analysis of the pion ratios from sect. 3.4, one 
favoring a very soft symmetry energy and the other a rather stiff one ra- Microscopic 
many-body results in the region up to 2po favor a behavior more similar to the flow 
experiment. 

I have attempted to give a brief overview of the determination of the nuclear 
symmetry energy in HIC’s. HIC’s are interpreted with transport theories and I have 
mentioned some of the challenges in such descriptions. Today a picture emerges where 
the information on the symmetry energy from HIC’s, nuclear structure, and neutron 
stars increasingly converges. But there are obviously also open problems, where a 
more thorough understanding of the mechanism and the analysis is needed. In the 
end it is neccessary to obtain a consistent picture for many observables in heavy ion 
collisions. 
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Figure 4: Summary of information on the density dependence of symmetry energy 
for below (left) and above (right) saturation density. On the left information is 
collected from nuclear structure and low energy HIC |4], while on the right the still 
not consistent information from higher energy collisions. The figure is discussed in 
the text. 
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